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Abstract 



We study the dynamics of a dissipative two-level, system driven by a 
monochromatic ac field, starting from the usual spin-boson Hamiltonian. The 
quantum Langevin equations for the spin variables are obtained. The ampli- 
tude of the coherent oscillations in the average position of the particle is stud- 
ied in the high temperature limit. The system exhibits quantum stochastic 
resonance in qualitative agreement with earlier numerical results. 
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I. INTRODUCTION 



Stochastic resonance is a nonlinear phenomena. It has been predicted and experimentally 
observed in damped classical systems[l-4]. It is characterized by a maximum in the response 
of the system to an external input signal as a function of noise strength at the input signal 
frequency. There have been some recent studies in stochastic resonance in damped quantum 
double- well systems as well [5-9] . Recently, Makarov and Makri have studied the phenomena 
in a dissipative two-level system numerically using an iterative path integral scheme and tried 
to understand the results analytically[8,9]. They demonstrate that it is possible to induce 
and maintain large amplitude coherent oscillations by exploiting the phenomena of stochastic 
resonance in quantum systems. In their treatment noise strength is varied, through the 
coupling to heat bath, to obtain stochastic resonance. The maximum in response is also 
obtained with respect to the driving field strength indicating breakdown of linear response 
theory. However, their analytical treatment involves several approximations. They start 
with the evolution of dynamical variables of a system using Heisenberg equation of motion 
for a given spin-boson (system-bath) Hamiltonian in the presence of a monochromatic field. 
In their weak coupling approximation the bath causes stochastic energy fluctuations which 
are unaffected by the dynamics of the two-level system . To obtain the steady state average 
position of the particle (a z ) they require steady state solution for the population difference 
between the two-level-system eigenstates (a x ) . The required expression for (a x ) is then 
taken from the solutions of the phenomenological optical Bloch equations in the rotating 
wave approximation. The expression for (cr z ) thus obtained with these approximations seems 
to be in good agreement with their numerical simulations. 

In the present work we study stochastic resonance analytically by systematically deriving 
the quantum Langevin equation for the system by eliminating bath variables [10- 13]. The 
response of the system to an oscillating field is obtained in the high temperature limit. This 
limit corresponds to treating the random force operator as a classical c-number variable. 
Our results are in qualitative agreement with earlier works. 
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In section II the derivation of quantum Langevin equation is presented. The spin Bloch 
equations are obtained and the stationary solution of the relevant variables is given in sec. 
III. The last section IV is devoted to results and discussion. 

II. The Quantum Langevin equation for system variables 

To obtain the Quantum Langevin equations for system variables we consider a symmetric 
two-level system interacting linearly with a bath of harmonic oscillators in the presence of 
a time dependent external monochromatic field Vocos(ut). The total Hamiltonian is given 
by [8,9] 

H = -hA a x + hu k ala k + J2 s 9k(ak + a\)a z + V cos(ujt)a z , (1) 

where a 's are the Pauli matrices. a k and a\ are annihilation and creation operators for 
bath variables and g k is the coupling constant and 2s is the distance between the two 
wells. The population difference between the two eigenstates of the two-level system, which 
are separated by 2hA , is given by (a x ). The average value of (a z ) represents the average 
position of the particle, (a z ) being +1 or -1 is equivalent to the particle in the right or the 
left well, respectively. 

The quantum equation of motion for any dynamical variable A can be obtained from the 
following evolution equation 

where [- , -] indicates commutator. The commutaton relations among the spin variables are 
given by 

[a x , a y ] = 2ia z , (3) 

and their cyclic permutations. Moreover, 

[Ojfc/,4] = 5 k y. (4) 

Using equations (l)-(4) one can readily write down the equations of motion for the 
variables as 
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and, 



da x 
~dt 

da y 



-2 E ^pKO) + al(t))<r v (t) ~ ^-cos(ut)a y (t), 



2A a z (t) + 2j2^ L (a k (t) + a\(t))a x (t) + ^cos(ut)a x {t) 



2V n 



d<? z 
dt 

dat 
~~dt 



U \ - S 09k /,x 

-iu k a k (t) - i——a z {t), 

h 



— = +iuj k a k {t) + i—^-a z {t). 
Equations (||) and (|9|) are linear and therefore can be explicitly integrated, to obtain 



a k (t) = a*(0)e-^ - i 



h Jo 



a 



al(t) =4(0)e^* + ^ / ^(tOe^^'^f. 

ft Jo 



(5) 
(6) 
(7) 
(8) 

(9) 

(10) 
(11) 



Where a k (0) and c4(0) are the bath operator values at the initial time t = 0. Substituting 

for a k (t) and a£(t) from equations (JTTJ) and ( pTTj) in equations and (0), we obtain the 

following Quantum Langevin equations for the system variables, 

da x 2s , s , s 2Vo , , , . 
— = — —F(t)a y {t) - — cos(ut)a y (t) 



1 

+ 2 



2«2 r 



ft 2 Uo 



sin(^-OK(0^'KOO 



+ E ^§rU(*) { / * sin M* - 0)^(0^} 



— f = 2A (j z (t) + F(t) + —cos{ujt)a x 

dt n Ti 



(12) 



£ ^P^(t) { /"* sin(c fc (t - i>*(0^ 



E ^ {jf MM* - t')a z (t')dt'(t) \ a x (t) 



(13) 
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Here F(t) is given by 

F(t) = J29k(a k (0)e-^ + 4(0)e^)- (14) 

k 

As the dynamical operators (afe(O) , a^(0)) of the bath are distributed in accordance 
with the statistical equilibrium distribution for given temperature T, F(t) is referred to 
as Langevin operator noise term [12,13]. The integrals in equations fljjD and (|i~3"D can be 
integrated by parts leading to 

+ i ? SS J C{^M + ^^}«w.-0)* 

-^E^K(0)^W + ^(t)a 2 (0)}co S (^t), (15) 

2 ft W(; 

= 2A o" 2 (t) + ^-F{t)a x {t) + ^cos(ut)a x (t) 
at n n 



+^ E "S^ R(0K(t) + <r«(fK(0)} <xw(w*i). (16) 

In equations (|i~5|) and (|i~6|) G{t — t') = Yj{^s^g\/h 2 uJk)cos{uJk(t — t')) represents a damping 
or memory kernel [10]. The last two terms in equation (|i~5D and ([16]) are transient terms 



which can be neglected if one assumes Ohmic spectral density for bath oscillators [10, 11], i.e., 

PH = \ E - "*) = owerpC--), (17) 

where a is a dimension less dissipation coefficient (or Kondo parameter). Notice that the 
form of the spectral density is not bounded from above and hence for physical reasons one 
introduces an upper cut-off frequency tu c , namely p(u) = auje~ < ^' UJc \ such that the frequency 
scale uo c is assumed to be much larger than the characteristic frequencies of the problem. 



With this form of spectral density one can readily show that the transient terms do survive 
up to a time scale (l/u c ), which can be made arbitrarily small and thus can be ignored [11]. 
Thus for the long time behaviour and for the Ohmic spectral density the equations (|1^) and 



(|T6|) get further simplified and we arrive at the equations of motion for spin variables as 



a j da z (t) f ^ da z {t) \ 

+ 2 -dT^ + ^-dT ' (18) 



^ = 2A a z {t) + ^F{t)a x {t) + ^cos(cot)a x 
at n n 

_£{^fl ff . (t) + , w *^l, (19) 



^ = -2A ot r„(f), (20) 

with the memory kernel G(t — t')=2a5(t — t'). 

Further simplifying equations (|I8|)-(|20|) and making use of the properties of spin opera- 
tors, we finally obtain the Langevin equations of motion for spin variables as 

da x f 2s 2V \ 2a A ( , 

— = {—Fit) - — cosM)j^(t) - — , (21) 

^ = 2A a z (t) + {^F(t) + ^cosM)}^(t), (22) 



^ = -2A a y (t). (23) 

The above Langevin equations involve the operator random force F(t). The statistical 
properties of F(t) can be obtained using the equilibrium distribution for bath variables, 

< G ^> = e ^]_ V ( 24 ) 
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where (3 = l/ksT. Using this and the Ohmic spectral density for the bath oscillators, the 
symmetrized autocorrelation for the operator valued random force F(t) is given by[12,13], 

^(F(t)F(t>) + F(t')F(t)} = Y.9l cosu k (t - t') coth(huj k /2k B T) (25) 

and the nonequal time commutator is given by 

[F(t),F(t')\ = -2iJ29 2 k smuj k (t-t')- (26) 
III. The Spin Bloch equation and their solutions. 

Owing to the operator nature of the random Langevin force F(t), it is difficult to solve 
for the expectation values of spin variables using equations - (p3|) . For simplification 
we make a first approximation in that the operator random force is treated as a classical 
c-number random variable. One can readily verify that in the classical limit [12,13], taking 
h — > 0, the nonequal time commutator of F(t) vanishes and the autocorrelation of the 
Gaussian random force F(t) becomes 

(F(t)F(t')) = V kT5(t-t'), (27) 

where i] is the friction coefficient and is related to Kondo parameter a through the following 
realation 

r) = (h/2s 2 )a. (28) 

The classical Markov approximation for F(t) is valid in the high temperature limit, 
which will become clear later. With the above approximation one can readily write down 
the equations of motion of spin variables averaged over the ensemble of realizations of random 
fluctuations F(t), whose autocorrelation is given by equation (p7|), with the help of Novikov's 
theorem [14-16]. We get 
d(a x (t)) 



dt 



S{a x (t)} - 2a A - A cos(ut)(a y (t)), (29) 



= 2A (a z (t)) - 5(a y (t)) + Acos(ujt)(a x (t)), (30) 
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= -2A (a,(t)>, (31) 

where 5 = 2aksT /h, and A = 2Vo/h. Equations fl2"9|) to ( plf ) represent the Bloch equations 
for spin variables. Unlike in the standard NMR situation, there is no relaxation term in the 
evolution of (cr z ) . This is because fluctuating environmental fields are exclusively along the 
z direction [10]. These equations are characterized by a single relaxation time r = <5 _1 = 
h/2akBT. In the absence of an external driving field, the system will relax asymptotically 
to equilibrium and the equilibrium value of population difference between the two levels 
separated by an energy value 2hA is given by (a x ) = (hAo/ksT) instead of tanh(/iA //i;_BT). 
This shows that our approximation of treating the operator random force as classical random 
c-number variable is valid for high temperature such that HAo/ksT <C 1. The stationary 
solution of Bloch equations in the presence of an external field can be found by using the 
method of harmonic balance. For this we assume stationary solutions for (cr z ) s and (cr x ) s to 
have the form 

{ <T z(t)) s = a cos (cut) + bsm(ujt) 

= (a z (t)} cos(ujt-<j)) (32) 

{cx{t)) s = V + ccos(atf) + dsin(ujt) (33) 

where a, b,y, c and d are constants to be determined. Substituting (|3~2"D and (|33"D in the Bloch 
equations ( ^9l) to (0) and using harmonic balance one readily gets an expression for the 
required steady state amplitude (cr z ) for the average position of the particle (a z ) as 



(a z } = Va^TV 

- AV A 2 V — 1 (qa\ 

11 K 2 - 4A 2 ) 2 + Aufa 2 + 2VM' 

Here in the above expression we have rescaled all the energy variables in terms of fc^T, i.e., 
V\ = Vo/ksT, Ai = HAa/hsT and uj\ = huj/ksT are dimensionless variables. Note that 
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in the absence of driving field the amplitude (cr z ) Q vanishes as expected. Away from the 
resonance {uj = 2A ), the oscillation amplitude for high frequency field scales as (1/uj 2 ). In 
the limit of low frequency (static) driving uj — > and for small Vq , i.e., in the adiabatic 
limit (cr z ) is independent of relaxation time r (= h/2akBT). These results are consistent 
with those obtained in ref. [8,9]. 
IV. Results and discussion 

In fig.(l) we have plotted the stationary amplitude of the average position of the particle 
(cr z ) as a function of the dissipation coefficient (or Kondo parameter) a, for various values of 
the dimensionless field amplitude V\ (= Vo/ksT). We have restricted to the case of resonant 
condition uj = 2A . The solid curve is for V\ = 0.01, long dashed curve for V 1 = 0.04 
and small dashed curve for V 1 = 0.07. We see that all these curves exhibit maxima at an 
optimum value of a, which depends on A^Vi and oj\. The maximum value of the peak 
(M p ), however, is independent of field strength and is given by M p = y/2(Al/ui), but the 
position of the peak (a m =y / 2l / 1 2 c<jf — (uj 2 — 4Af) 2 /2c^i) shifts towards higher values of a as 
we increase V\. The independence of the value of peak maxima on the field strength is also 
noted in reference [8,9], for resonant case. Our expression for M p is valid for any frequency. 
The occurrence of the peak or the maximum in the (<J Z ) as a function of coupling strength 
a is attributed to stochastic resonance in quantum two-level systems. This is a result of 
cooperative phenomena between various competing mechanisms of energy exchange between 
the two-state system with thermal bath and the external driving field. 

In fig. (2) and fig. (3) we have plotted (cr z ) versus UJ\{= hu/ksT). For this we have 
taken Vi(= V /k B T)=QM and A x (= A/A; B T)=0.2. For fig. (2) a = 0.05 and for fig.(3) 
a = 0.01. Fig. (2) and (3) should be compared with fig. (7a) and fig. (7b) of respectively 
ref. [9] for their qualitative behaviour. In fig. (2) (cr z ) exhibits a maximum at the resonant 
frequency uj = 2A , whereas two maxima are seen in fig. (3) at off resonance frequency values. 
This indicates stochastic resonance can be obtained even for the off resonance conditions. 
This further indicates that the stochastic resonance is indeed a bonafied resonance [17,18]. 

In fig(4) we have plotted (cr z ) as a function of external driving field amplitude V\{= 
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Vo/ksT) for various values of the dissipation parameter, the frequency of the field is lo\ = 
2A 1 = 0.4. The solid, long dashed and small dashed curves are for a = (0.02), (0.04) 
and (0.06), respectively. We notice that the response of the system, i.e., (o" 2 ) , initially 
increases with the field amplitude and attains a maximum value at a particular value of V\ 
depending on a and on other parameters. After exhibiting the maximum (cr z ) decreases as 
we increase V\ further. The maximum value (M p ) of the peak in the (cr z ) for a given value 
of Ai, loi and a occurs at the field amplitude value V\ = y/ iK- 4 ^^" 1 " 4 " anc i j s equal 

— A 2 

to M p = v2(^) which is independent of a, the Kondo parameter. The linear response 
theory is valid for small V\ to the left side regime of the maximum. In these plots stochastic 
resonance manifests as a breakdown of linear response theory, thus bringing out the non 
linear nature of the problem explicitly. 

In fig. (5) we have plotted (<J Z ) versus A x (= hAo/ksT), for given uj\{= hu/k B T)=0A 
and a small a=0.04. The maximum appears well within the range of acceptable values of 
Ai < 1, i.e., hA <C k B T. For very small values of a we do get maximum response at two 
values of Ai(^ 2ui) indicating stochastic resonance at off resonance condition (similar to 
the observation made in fig. (3)). 

In conclusion we have derived quantum Langevin equation for a dissipative two-level sys- 
tem, driven by monochromatic ac field, starting from microscopic spin-boson Hamiltonian. 
The equations of motion for the average values of the spin variables are obtained in the high 
temperature limit. We have obtained an analytical expression for the amplitude of coherent 
oscillation in the average position of the particle, which exhibits stochastic resonance with 
respect to various parameters in the problem. These results are in qualitative agreement 
with the results obtained by Makarov and Makri using the numerical method of iterative 
path integral scheme [8,9]. 
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Figure caption 

Fig.l Plot of the stationary amplitude of average position (a z ) versus disipation co- 
efficient (Kondo parameter) a for a fixed value of cj=2A =0.4. For this case V\{= 
Vo//cb7 1 )=0.01 (solid curve), 0.04(long dashed curve) and 0.07(small dashed curve). 

Fig2. Plot of the stationary amplitude of average position (cr z ) versus oj\ (= hu /kBT) 
for Vi=0.04, Ai=0.2 and a=0.05. 

Fig3. Plot of the stationary amplitude of average position (a z ) versus io\ (= fiuo/hsT) 
for 14=0.04, Ai=0.2 and a=0.01. 

Fig. 4 Plot of the stationary amplitude of average position (cr z ) Q versus V±(= Vo/ksT) for 
a fixed value of oj=2A =0.4. For this case a=0.02 (solid curve), 0.04(long dashed curve) 
and 0.06(small dashed curve). 

Fig. 5 Plot of stationary amplitude of average position (cr z ) Q versus A 1 = (hA /k B T) for 
a fixed value of Vi=0.04, ui =0.4, and a =0.04. 
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